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Contract DA-49~186-502-ORD-860 
Revised Scope of Work - ORDTL-CB-762 

HARRY. DIAMOND LABORATORIES 
Task 1 - Energy Transformation Systems 

Conduct theoretical and experimental investigations to define 
the mechanism and control parameters of gas diffusion electrodes to 
fuel cells. Work will include, but not be limited to 

1. Establish a mathematical model of the overall diffusion 
process including terms for contributions of known mechanisms in 
porous media. 

2. Develope expressions for contributing mechanisms which 
contain the physical properties of the reactants and products, and 
the physical properties and geometrical parameters of the electrode. 

3. Conduct experimental work to verify the mathematical model. 


1.0 SUMMARY 


A mathematical model is developed for the mass transport of 

a reacting species through a liquid to an electrode surface. The 

resistance is considered to be concentrated in a stagnant film of 

thickness 6. By considering no bulk flow of the non-reacting 

species the following equation is derived. 

. -eng 
^m 6 

If bulk flow of component B is, considered, equation (2) is 
obtained. 


( 1 ) 


m °AB ,, m+1 o. 


( 2 ) 


In equations (1) and (2) the symbols have the following significance: 
i^ = current density supportable by mass transport 
c = concentration of solvent - •- 
n = number of g,. equivalent per gram mole 
F = Faraday 

= diffusivity of A in B 

= concentration of A in bulk * * 

A 

6 = film thickness 

m = moles of B produced per mole of A 

Calculations are made for the case of oxygen (A) in water (B) 

under various concentrations and temperatures. For a film thickness 

of 10 ^ cm, a temperature of 300^K, and a partial pressure of oxygen 

-2 2 

of 1 atm, i^ is equal to 1.03 x 10 amps/cm . 

For the case examined there is no significant difference between 
equations (1) and (2), 


1 





2 


It is found that current densities are highest for low temper¬ 
atures and high partial pressures of At constant pressure of 0 ^, 
the current density exhibits a minimuni at approximately 80^0, 




2,0 INTRODUCTION 


One possible method of obtaining electrical energy directly 
from chemical energy is in a fuel cell with submerged electrodes. 
In this type of cell two compartments, the anode and the cathode 
compartments, might be filled with aqueous solutions saturated 
with an oxidizable gas and oxygen respectively. The gas may be 
bubbled through the solution in a compartment to keep it satur¬ 
ated and well mixed. In this type of cell the rate-controlling 
step could be the rate of diffusion of the reactants to the 
electrode surface. 

The purpose of this report is to calculate the current 
density sup^portable by mass transfer to such submerged electrodes. 



3.0 MATHEMATICAL MODEL 


Since the processes in each electrode compartment are similar it is only 
necessary to develope the mathematical expressions for one of them, A solution 
for the other compartment results when the proper values of the several variables 
are inserted in the equations. 

To develop a model in general terms, consider A diffusing through B. 

Since the liquid, B, is assumed to be well mixed, the concentration of A at 

any point in the bulk liquid may be assumed to be constant . We will also 

assume that all A which reaches the electrode surface will react instantaneously, 

so that the concentration of A at the electrode surface is zero. Because the 

concentration difference between the bulk liquid and the electrode surface is 

the only driving force and the electrode may be thought of as having infinite 

area as compared to the thickness of liquid through which reactant molecules must 

diffuse, we are concerned only with one-dimensional transport. This concentration 

difference may be assumed to exist over a small distance, 6, which may be thought 

of as a stagnant film. One may derive an expression for the flux, flow per unit 

time per unit area, of A arriving at the electrode surface by diffusion through 

this stagnant film. A schematic diagram is shown in Figure 1. 
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Figure 1 

Schematic Diagram of Cathode Compartment 
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Writing a shell mass balance (1) 

dx. 


"a ■ -'V 3r + *a<“a ♦ ’*B> 


3 . 0-1 


In equation ’3.0-1 the first term on the right hand side is the flux 

resulting from diffusion and the second term is the flux resulting from bulk 

flow. Before this equation can be solved a relation between N. and N_ must 

A D 

be known. 


3.1 No Bulk Flow of Component B 


If we assume N_ to be equal to zero, equation 3.0-1 now becomes 
jd 

- V = -^°AB IT 


3.1.1 


Two solutions to this equation are possible. If we assume x. « 1 and 


D._ not a function of concentration we obtain 
AB 


r ■ r- 




3.1-2 


N . .^A 

\ 6 


3.1-3 


Equation 3.1-3 is the Pick’s law equation (2) for steady state diffusion 
of A through stagnant B. 


If X. = 1 then the solution to equation 3.1-1 is 
* 




3.1-4 


n 1 1 

N, = —;— In 


A 6 o 

^ - ^A 


3.1-5 
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3.2 Bulk Flow of Component B 

If we consider the bulk flow of B, we take into account the removal 
of the reaction product, assuming that A reacts to form B, as it would 
in the case of oxygen and water. 

A = mB. 3.2-1 

N = + m 3.2-2 

A B 


Equation 3.0-1 now becomes 

m + 1 




1 - 


m A 


= -cD 


AB dz 


3.2-3 


and integrating with not a function of concentration 

„ m , ,, m + 1 o\ 

~ —ITT "x— (1 - - x .*) 

A. m + 1 0 m A 


3.2-4 


3.3 Numerical Calculations 

i^, the current density supportable by mass transport is directly 
proportional to the flux. 


i = nFN. 3.3-1 

m A 

4 

F = 9.65 X 10 amp sec/gram equivalent 
n = number of gram equivalents/gram mole 

For investigational purposes we will assume that A is oxygen and 
B is water. 

Applying Henry’s law to this situation 



3.3-2 


we obtain for equation 3.1-3 in terms of current, density 





7 


nF cD- „ ^ P. 

O2 ~ ^2° ^2 


3.3-3 


The c in equation 3.0-1 is the concentration of water (B) in the 
solution of water and oxygen (A)- However since this is a very dilute 
solution of oxygen, it is taken equal to the concentration of water in 
water, or 0.0556 moles/cc. 

For low concentrations of A in B, Bird (1) reports that the 
diffusivity is directly related to the absolute temperature. 


n ^ tP — 


3.3-4 


Brown (3) reports the value of ^ at 20°C to be 2,08 x 10 ^ 

2 

cm /sec. With 4 gram equivalents of oxygen per mole equation 3.3-3 now 


becomes 


-3 T ’ ^2 

i = 1.52 X 10 7 

m 0 . H 


3.3-5 


There is very little information on the values of 6 in the liter¬ 
ature. To show the relationship of i with temperature and pressure, let 

m 

us introduce a new variable O = (i )( 6 ). Calculations of $ are then 

m 

made at various temperatures and pressures. H, the Henry's law constant, 
is a function of temperature and pressure. Values from Perry's Chemical 
Engineers Handbook s page 675, Tables 21, and 22 may be used.. 

A normalizing factor, at 300^K and 1 atmosphere partial pressure 
of oxygen is used to obtain a normalized current density, which values 
are shown in Table 1, and Figures 2 and 3. This value of is 1.03 x 
10 ^ cm amps/cm^. 
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Table 1 

NORMALIZED CURRENT DENSITY AT A SUBMERGED ELECTRODE 


t?C 

,atm 

T,K X 

0.1 

0.5 

I.O 

6.0 

10.0 

7 

280 

0.138 

0.697 

1.38 

6.70 

10.07 

17 

290 

0.116 

0.583 

1.161 

5.874 

9.462 

27 

300 

0.1 

0.502 

1 

5.202 

8.475 

37 

310 

0.089 

0.446 

0.892 

4.776 

•7.780 

47 

320 

0.083 

0.415 

0.830 

4.529 

7.354 

57 

330 

0.080 

0.398 

0.798 

4.350 

7.085 

67 

340 

0.078 

0.388 

0.776 

4.260 

6,951 

77 

350 

0.076 

0.390 

0,761 

4.22 

6.906 

87 

360' 

0.076 

0.381 

0,762 

4.215 

6.951 

100 

373 

0.0789 

0.394 

0.789 

4.318 

7.085 


Some discussion should be made of 6, the effective film thickness. 
This is perhaps best thought of as a hypothetical layer whose thickness 
is such that it gives the proper mass transport, as observed experimen¬ 
tally, given a fixed concentration gradient. This film is simply that 
region where the concentration gradient is at a maximum. It is a cal¬ 
culated quantity, not easily measured experimentally, and therefore, few 

-2 

values exist in the literature. It is reported that 2.5 x 10 cm(4) 

is not an unreasonable figure for this film* Other references (5,6) 

-3 -2 

gives values of 1 x 10 cm for stirred solution and 5 x 10 cm for 


unstirred solutions 
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Calculating current density at 1 atmosphere oxygen partial pressure 
o "3 

and 300 K,, with a 6 of 10 cm thick, we obtain a value of 1.03 x 10 

amps/cm^. Other values can be obtained, for any 6, by using ^ and the 
figures in Table 1. The thinner this film can be made the higher the 
current density possible. 


Current densities are highest for low temperatures (7 C) and high 
partial pressures of 0^ because the solubility of 0^ is greater under 
these, cirexamstances. At constant pressure, the current density exhibits 
a minimum at approximately 80^C and then shows a slight increase due to 
the decrease of 0^ solubility and the increase of 0^ diffusivity with 
increasing temperature. 

If we consider the bulk flow of water (B) m, in equation 3.2-2 is 


equal to -2 and equation 3.2-4 becomes 


1 = +2 D ln(l - -r x.) 

m 0 ^2~^2^ 2 A 


3.3-5 


For very small y, ln(l + y) is approximately equal to y, and 
substituting equation 3.3-2 and 3.3-4 

P 

- = nEc yP 1 _^ 

6 ^O^-H 0 TO H 


3.3-6 


i = 1.52 X 10 ^ 3.3-7 

m o H 

By considering bulk flow of the solvent as in equation 3.2-4, we find 
that for the case of 0^ in H^O the current density is uneffected. An 
exact calculation without the approximation ln(l + y) = y would show 
a decrease in current density for the case of bulk flow of the solvent. 







SYMBOLS 


constant concentration of solvent 
diffusivity of o-xygen in, water 

Henry*s law constant atms/mole 

current density supportable by mass transport 

4 . 

conversion factor 9,65 x 10 ampere-seconds/isquivalent 

2 

mass flux of ^2 ^oles/sec cm 

partial pressure of oxygen in atm 

absolute temperature, 
mole fraction, 
thickness of film, cm 

2 

current density times film thickness ( 6 )(i^)(amps/cm )(cm) 
number of g. equivalent per gram mole 
moles of B produced per mole of A 
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